Neutrino parameters from matter effects in Pee at long baselines 
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We show that the earth matter efTects in the Uc Vc survival probabihty can be used to cleanly 
determine the third leptonic mixing angle S13 and the sign of the atmospheric neutrino mass squared 
diflterence, Am^i, using a /3-beam as a u^. source. 
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Determination of the third neutrino mixing angle ^13, 



the sign of Am: 



the three CP phases, and 



the absolute neutrino mass scale are necessary for recon- 
struction of the neutrino mass matrix, which will have 
important consequences for nuclear and particle physics, 
astrophysics and cosmology. The Ve — > transition 
probability Pe^i , has been identified as the "golden chan- 
nel" [l[ for measuring the Dirac phase 5c sgn{lS.m'i^^) 
and in long baseline accelerator based experiments. 
However, this strength of the golden channel also brings 
in the well-known problem of parameter "degeneracies" , 
where onegets multiple fake solutions in addition to the 
true one Various ways to combat this vexing issue 
have been suggested in the literature, including combin- 
ing the golden channel with the "silver" (Per) H and 
"platinum" (P^e) channels. While each of them would 
have fake solutions, their combination helps in beating 
the degeneracies since each channel depends differently 
on 5cp, sg7i(Am|i) and ^13. In this letter, we propose 
using the u,, — > Vf, survival channel, Pee^ which is inde- 
pendent of 5c p and the mixing angle ^23- It is therefore 
completely absolved of degeneracies and hence provides a 
clean laboratory for the measurement of sgn{Am\i) and 
9 13. This gives it an edge over the conversion channels, 
which are infested with degenerate solutions. 

The Pee survival channel has been extensively consid- 
ered for measuring ^13 with Ve produced in nuclear reac- 
tors [3] and with detectors placed at a distance ~ 1 km. 
Reducing systematic uncertainties to the sub-percent 
level is a prerequisite for this program and enormous 
R&D is underway for this extremely challenging job. For 
accelerator based experiments, the survival channel, Pee^ 
has been discussed with sub-GeV neutrinos from a (3- 
beam source at CERN and a megaton water detector in 
Prejus at a basehne of 130 km jS, 1 .'i . However, no sig- 
nificant improvement on the 6'i3 limit was found in [5] 
for a systematic error of ^ 5%. This stems mainly 
from the fact that in these experiments one is trying to 
differentiate between two scenarios, both of which pre- 
dict a large number of events, differing from each other 
by a small number due to the small value of 6113. Also, 
since sgn{/!^m\i), is ascertained using earth matter ef- 
fects, there is no hierarchy sensitivity in these survival 



channel experiments due to the the short baselines in- 
volved. 

In this letter, we emphasize on the existence of large 
matter effects in the survival channel. Pee, for an exper- 
iment with a very long baseline. Recalling that Pee = 
1 — Pef_i — Per and since for a given sgn{Am^i) both Pe^ 
and Per will either increase or decrease in matter, the 
change in Pee is almost twice that in either of these chan- 
nels. Using the multi-GeV Ve flux from a /?-beam source, 
we show that this large matter effect allows for signifi- 
cant, even maximal, deviation of Pee from unity. This, 
can thus be a convenient tool to explore ^13. This is in 
contrast to the^ reactor option or the /3-beam experimen- 
tal set-up in [5:], where increasing the neutrino flux and 
reducing the systematic uncertainties are the only ways 
of getting any improvement on the current ^13 limit. We 
further show, for the first time, that very good sensitivity 
to the neutrino mass ordering can also be achieved in the 
Pee survival channel owing to the large matter effects. We 
discuss plausible experimental set-ups with the survival 
channel and show how the large matter effect propels this 
channel, transforming it into a very useful tool to probe 
S(7n(Am3j^) and ^13 even with relatively large room for 
systematic uncertainties. 

For simplicity, we start with one mass scale domi- 
nance (OMSD) and the constant density approximation. 
OMSD implies setting Am^i = 0. Under these condi- 
tions, for neutrinos of energy E and traveling through a 
distance L, 
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[l.27(Am2j™L/P] , (1) 



where the mass squared difference and mixing angle in 
matter are respectively 



(Am^i) 
sin 29 



13 



( Am^i cos 26'i3 - 
sin 26*13 Am^i/(ATO^i)'" 



(Am§i sin 26*13)2 
(2) 



and A = 2^/2G pneE originates from the matter poten- 
tial Here, Ue is the ambient electron density. From 
Eq. ([1]), the largest deviation of Pee from unity is 
obtained when the conditions (i) sin^ 26*5^ — 1 and (ii) 
sin^ [l.27(Am§i)™L/i?] = 1 are satisfied simultaneously. 
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The first condition is achieved at resonance which is ob- 
tained for A = Ato3]^ cos 26*13. This defines the resonance 
energy as, Eres — ^^31 cos 20i3/2-\/2Gi?ne. The second 
condition gives the energy where the (Am|i)™ driven 
oscillatory term is maximal, 

_ 1.27{Aml,rL 
E^a.- (2p+ 1)^/2 'P-0,l>2.. (3) 

Maximum matter effect is obtained when Eres = -^'max 



0, which gives. 



(2p+ l)7r 5.18 X 10^ 
tan 29^ 3 



km gm/cc. (4) 



p is the matter density in gm/cc. This is the distance 
where Pee — 0. Although both Eres and E^^^ depend 
on the value of Am|j^, the distance at which we get the 
maximum matter effect is independent of Amg^. How- 
ever, it is controlled by ^13 very sensitively. For aver- 
age earth matter densities obtained using the Prelimi- 
nary Reference Earth Model (PREM), one can find the 
typical distances at which the above conditions are sat- 
isfied for various values of sin^ 2^13 [Sj] . For instance for 
p = and sin^ 26'i3 = 0.2 and 0.1, these distances are 
7600 km and 10200 km respectively. For higher values of 
p the distance exceeds the earth's diameter for 9iz in the 
current allowed range. Using {pL)"^°-^ corresponding to 
the PREM profile, from Eq. @ one can estimate that 
the condition of maximal matter effects inside the earth's 
mantle is satisfied only for sin^ 2^13 ^ 0.09. 

Under OMSD, the matter conversion probabilites 



Pe 



^23 sin^ 26*5^ sin^ [l.27(Amii)"L/£;] , where 



Y23 — sin 6*23 for x — p and 1^23 — cos ^23 for x — t. 
The maximum matter effect condition in the conversion 
channels is also given by Eq. ([4]). However, there are 
suppression factors, sin^ ^23 for Pe/x and cos^ ^23 for Per , 
not present in Pee- Moreover, since Pee does not contain 
023, the octant ambiguity as well as parameter correla- 
tions due to uncertainty in ^23 are absent. In addition, 
as mentioned earlier, the Pee channel does not contain 
the CP phase, 5c p- Both of these remain true in the 
presence of non-zero Am|i In our numerical work, 
we solve the full three flavour neutrino propagation equa- 
tion assuming the PREM [13] profile and keep Am|i and 
sin^ 9x2 fixed at their present best-fit values of 8.0 x 10"^ 



eV and 0.31 respectively [ll|. We assume the true value 



of |Am§i 



2.5 X 10-3 eV^. 
In Fig. [1] we plot Pee as a function of energy, at four 
different L and for three values of sin^20i3. The plots 
confirm that maximal matter effects come at i ~ 10000 
km and L ~ 7500 km for sin^ 2^13 — 0.1 and 0.17 respec- 
tively for the normal hierarchy (NH). For the inverted 
hierarchy (IH) there is no significant matter effect for 
Ve- This large difference in the probabilities for NH and 
IH can be exploited for the determination of sgn{l\m'^^) . 
Further, since the matter effect is a sensitive function of 
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FIG. 1: Pee in matter plotted versus neutrino energy. Thick 
(thin) lines are for normal (inverted) hierarchy. 
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FIG. 2: Events in 5 years vs. sin^ 2^13 for the normal (dashed 
line) and inverted (solid line) hierarchy for L = 7500 km and 
7 = 500. The inset shows the same but for 7 = 250. 



013 it may also be possible to obtain information on this 
angle. We can also see that for a given value of sin^ 2^13 
( ^ 0.09) and P, the matter effect increases (almost lin- 
early) with L, until the L for maximal matter effect is 
reached, beyond which matter effect falls. For values of 
sin^ 2013 < 0-09 the condition for maximum matter ef- 
fect is not met inside the earth's mantle and hence the 
matter effect and sensitivity to both hierarchy as well as 
013 increase with L. 

In what follows, we will show how, in a plausible exper- 
iment, one can use this near-resonant matter effect in the 
survival channel. Pee, to constrain 0i3 and sgn(l\m\i). 
Fig. [Ushows that the requirements for such a program in- 
clude a Ve beam, baselines of at least a few thousand km 
and average energies around 6 GeV. The detector should 
be able to observe e~ unambiguously at these energies. 

Pure Vejve beams can be produced from completely 
ionized radioactive ions accelerated to high energy decay- 
ing through the beta process in a storage ring, popularly 
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known as /3-beams 12|, |l3|, [ij] . The ions considered as 
possible sources for beta beams are ^^A'^e and for 
and ^He and ^Li for Pg. The end point energies of ^He 
and ^^iVe are 3.5 MeV while for and ^Lz this can 
be larger ~ 13-14 MeV [ISl]. For the Lorentz boost factor 
7 — 250(500) the ^He and ^^Ne sources have peak energy 
around ~ 1(2) GeV whereas for the ^B and ^Li sources 
the peak occurs at a higher value around ~ 4(7) GeV. 
Since the latter is in the ball-park of the energy neces- 
sary for near-resonant matter effects as discussed above, 
we will work with ^B i^Li) as the source ion for the v^. 
(Pe) /3-beam and 7 = 250 and 500. A CERN /3-beam 
facility with 7 ~ 250 should be possible with the exist- 
ing SPS, while 7 < 500 could be achieved with upgrades 
of the existing accelerators. The Tevatron is also being 
projected as a plausible accelerator for the /3-beam. 

Water Cerenkov detectors have excellent capability of 
separating electron from muon events. Since this tech- 
nology is very well known, megaton water detectors are 
considered to be ideal for observing /?-beams. Such de- 
tectors do not have any charge identification capacity. 
But in a /3-beam, the (3~ and emitters can be stacked 
in different bunches and the timing information at the 
detector can help to identify the e~ and e"*" events [l6l |. 

It is well known that there are no beam induced back- 
grounds for /3-bcams. In this experimental set-up, the 
process u,, ^ i't ^ ^ could mimic the signal. We 
have checked that the background to signal ratio for these 
events in the relevant energy range is ~ 10~^ and can be 
neglected for the disappearance mode. e~ events from 
K and tt~ decays are also negligible. The atmospheric 
background can be estimated in the beam off mode and 
reduced through directional, timing, and energy cuts. 

Proposals for megaton water detectors include UNO 
[17] in USA, HyperKamiokande in Japan and MEM- 
PHYS J^] in Europe. If the /3-beam is produced at 
CERN, then basehnes in the range 7000-8600 km would 
be possible at any of the proposed locations for the 
UNO detector. Likewise, if the /9-beam source be at 
FNAL, then the far detector MEMPHYS would allow 
for L = 7313 km. HyperKamiokande could also be con- 
sidered as the far detector and in that case L = 10184 
(9647) km if the source be at FNAL (CERN). 

For our numerical analysis we use the standard 
technique with xlotai = xluU < ^prtor, 
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sider 2% /3-beam flux normalisation error and 2% error 
for detector systematics which is more conservative than 
the value of 2% usually considered in current literature 
For the full definition of y?„;; and details of 
our statistical analysis, we refer to [20|. The prospective 
"data" is generated at the "true" values of oscillation 
parameters, assuming 440 kton of fiducial volume for the 
detector with 90% detector efficiency, threshold of 4 GeV 
and energy smearing of width 15%. For the i>e /3-beam we 
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FIG. 3: Sensitivity to hierarchy for L — 7500 (solid line) and 
10000 km (dashed line) and 7 = 250, as a function of true 
value of sin^ 2613 . 



0.2 
0.1 
0.05 

0.02 

0.01 
0.006 
0.1 
0.05 

0.02 
0.01 
0.004 



1 














H 1 1 


■ \ — 

7 = 


\ \ 

500(30 


H 

y 


= 250(30) 




500(90%) 


y 


= 250(90%) ■ ' 



























2000 3000 4000 5000 6000 7000 8000 9000 10000 
L (km) 

FIG. 4: The upper panel shows the smallest value of true 
sin^ 2013 at which normal and inverted hierarchy can be 
distinguished while the lower panel gives the sensitivity to 
sin"^ 26i3 at various baselines at 90% and 3cr C.L., for two 
values of 7. 



have assumed 1.1 x 10^^ useful ^B decays per year and 
show results for 5 years of running of this beam. The 
number of events as a function of sin^ 26*13 at L = 7500 
km with a 7 = 500 i^e /3-beam is shown in Fig. [5] for NH 
and IH. The inset in Fig. [2] shows the number of events 
in 5 years expected from a lower 7 = 250. We have used 
the neutrino- nucleon interaction cross-sections from [2H . 

In Fig. [3] we show the sensitivity (ncr, n = \f^) of 
the survival channel to the neutrino mass ordering for 
L = 7500 and 10000 km and 7 = 250. If the true value 
of sin^ 2013 = 0.05, then one can rule out the inverted 
hierarchy at the 4.8cr (5. Oct) level with L = 7500 (10000) 
km. For L = 7500 (10000) km, the wrong inverted hier- 
archy can be disfavored at the 90% C.L. if the true value 
of sin^ 26*13 > 0.03 (0.025). The sensitivity improves sig- 
nificantly if we use 7 = 500 instead of 250, since (i) the 
flux at the detector increases, and (ii) the flux peaks at 
E closer to 6 GeV, where we expect largest matter effects 
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[20|. For 7 = 500, the inverted hierarchy can be disfa- 
vored at 2.6a{3.8a) for a lower value of sin^ 29i3 = 0.015. 
Minimum values of sin^ 2^13 at which the inverted hier- 
archy can be ruled out at 90% and 3a C.L. for different 
values of L are shown in the upper panel of Fig. 2] for 
7 — 250 and 500. From the figure one can see that for 
7 = 500 and L = 7500 (10000) km the wrong inverted 
hierarchy can be disfavored at the 90% C.L. if the true 
value of sin^ 26*13 > 1.0 x 10"^ (8.0 x 10"^)^ If instead we 
use a total systematic error of 5% then we get the above 
sensitivity limits as sin^ 26*13 > 1-6 x 10"^ (1.2 x 10"^) 
at 90% C.L. 

If the true value of sin^ 26*13 turns out to be smaller 
than the sensitivity reach shown in the upper panel of 
Fig. [4] for a given L, then it would not be possible to 
determine the hierarchy at the given C.L. However, we 
would still be able to put better constraints on sin^ 26*13 
itself. The lower panel of Fig. [4] demonstrates as a func- 
tion of L the sensitivity to ^13, i.e., the minimum value 
of sin^ 2^13 which can be statistically distinguished from 
sin^ 2013 = at 90% and 3a C.L. Both Figs. [3] and 0] 
show that the sensitivity improves with L, even though 
the flux falls as This results from matter effects in- 

creasing with L, as noted before. For L = 7500 (10000) 
km, we can constrain sin^ 2^13 < 6.3 x 10~^ (4.3 x 10~^) 
at the 90% C.L. for 7 = 500. For a 5% systematic error 
the above numbers are changed to sin^ 2^13 < 1.0 x 10^^ 
(7.3 X 10-3). 

How does this compare with alternate possibilites ? If 
the energy can be reconstructed accurately, then the re- 
sult can be improved further. For instance, for L=7500 
km, if one could preferentially select the energy in the 
range 5.0-7.5 GeV, then the normal and inverted hierar- 
chy would be differentiated for sin^ 2^13 — 7.47 x 10^^ at 
90% C.L. for 7 = 500. 

Using the Pep, channel and a /3-beam source at a dis- 
tance of 7152 km (CERN-INO) from the detector, the NH 
can be distinguished from IH at 90% C.L. for 7 = 500 
if sin^ 26*13 > 7.7 x 10"^ [10]. For a neutrino factory at 
7500 km, the wrong IH can be discarded at the 3a C.L. 
if sin^ 26*13 > 3 X 10-4 

We have presented our results using a v,, beam and 
assuming NH to be the true hieararchy. Similar results 
can also be obtained with a beam for IH. It is also 
possible to run both beams simultaneously. 

In conclusion, we propound the possibility of using 
large matter effects in the survival channel, P^e, at long 
baselines for determination of the neutrino mass ordering 
(sgrt(Am|]^)) and the yet unknown Icptonic mixing an- 
gle 6*13. Matter effects in the transition probabilites P^^ 
and Per act in consonance to give an almost two-fold ef- 
fect in the survival channel. In addition, the problem 
of spurious solutions due to the leptonic CP phase and 



the atmospheric mixing angle 823 does not crop up. The 
development of /3-beams as sources of pure Ue/ve beams 
enables one to exclusively study the Pee survival prob- 
ability and adds a new direction to the prospects of a 
future /3-beam. 
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